3.2.6 \(\int \frac {A+B x^2}{x^3 (a+b x^2+c x^4)} \, dx\) [106]

Optimal. Leaf size=112 \[ -\frac {A}{2 a x^2}-\frac {\left (A b^2-a b B-2 a A c\right ) \tanh ^{-1}\left (\frac {b+2 c x^2}{\sqrt {b^2-4 a c}}\right )}{2 a^2 \sqrt {b^2-4 a c}}-\frac {(A b-a B) \log (x)}{a^2}+\frac {(A b-a B) \log \left (a+b x^2+c x^4\right )}{4 a^2} \]

[Out]

-1/2*A/a/x^2-(A*b-B*a)*ln(x)/a^2+1/4*(A*b-B*a)*ln(c*x^4+b*x^2+a)/a^2-1/2*(-2*A*a*c+A*b^2-B*a*b)*arctanh((2*c*x
^2+b)/(-4*a*c+b^2)^(1/2))/a^2/(-4*a*c+b^2)^(1/2)

________________________________________________________________________________________

Rubi [A]
time = 0.16, antiderivative size = 112, normalized size of antiderivative = 1.00, number of steps used = 7, number of rules used = 6, integrand size = 25, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.240, Rules used = {1265, 814, 648, 632, 212, 642} \begin {gather*} -\frac {\left (-2 a A c-a b B+A b^2\right ) \tanh ^{-1}\left (\frac {b+2 c x^2}{\sqrt {b^2-4 a c}}\right )}{2 a^2 \sqrt {b^2-4 a c}}+\frac {(A b-a B) \log \left (a+b x^2+c x^4\right )}{4 a^2}-\frac {\log (x) (A b-a B)}{a^2}-\frac {A}{2 a x^2} \end {gather*}

Antiderivative was successfully verified.

[In]

Int[(A + B*x^2)/(x^3*(a + b*x^2 + c*x^4)),x]

[Out]

-1/2*A/(a*x^2) - ((A*b^2 - a*b*B - 2*a*A*c)*ArcTanh[(b + 2*c*x^2)/Sqrt[b^2 - 4*a*c]])/(2*a^2*Sqrt[b^2 - 4*a*c]
) - ((A*b - a*B)*Log[x])/a^2 + ((A*b - a*B)*Log[a + b*x^2 + c*x^4])/(4*a^2)

Rule 212

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))*ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x]
/; FreeQ[{a, b}, x] && NegQ[a/b] && (GtQ[a, 0] || LtQ[b, 0])

Rule 632

Int[((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(-1), x_Symbol] :> Dist[-2, Subst[Int[1/Simp[b^2 - 4*a*c - x^2, x], x]
, x, b + 2*c*x], x] /; FreeQ[{a, b, c}, x] && NeQ[b^2 - 4*a*c, 0]

Rule 642

Int[((d_) + (e_.)*(x_))/((a_.) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> Simp[d*(Log[RemoveContent[a + b*x +
c*x^2, x]]/b), x] /; FreeQ[{a, b, c, d, e}, x] && EqQ[2*c*d - b*e, 0]

Rule 648

Int[((d_.) + (e_.)*(x_))/((a_) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> Dist[(2*c*d - b*e)/(2*c), Int[1/(a +
b*x + c*x^2), x], x] + Dist[e/(2*c), Int[(b + 2*c*x)/(a + b*x + c*x^2), x], x] /; FreeQ[{a, b, c, d, e}, x] &
& NeQ[2*c*d - b*e, 0] && NeQ[b^2 - 4*a*c, 0] &&  !NiceSqrtQ[b^2 - 4*a*c]

Rule 814

Int[(((d_.) + (e_.)*(x_))^(m_)*((f_.) + (g_.)*(x_)))/((a_.) + (b_.)*(x_) + (c_.)*(x_)^2), x_Symbol] :> Int[Exp
andIntegrand[(d + e*x)^m*((f + g*x)/(a + b*x + c*x^2)), x], x] /; FreeQ[{a, b, c, d, e, f, g}, x] && NeQ[b^2 -
4*a*c, 0] && NeQ[c*d^2 - b*d*e + a*e^2, 0] && IntegerQ[m]

Rule 1265

Int[(x_)^(m_.)*((d_) + (e_.)*(x_)^2)^(q_.)*((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4)^(p_.), x_Symbol] :> Dist[1/2,
Subst[Int[x^((m - 1)/2)*(d + e*x)^q*(a + b*x + c*x^2)^p, x], x, x^2], x] /; FreeQ[{a, b, c, d, e, p, q}, x] &&
IntegerQ[(m - 1)/2]

Rubi steps

\begin {align*} \int \frac {A+B x^2}{x^3 \left (a+b x^2+c x^4\right )} \, dx &=\frac {1}{2} \text {Subst}\left (\int \frac {A+B x}{x^2 \left (a+b x+c x^2\right )} \, dx,x,x^2\right )\\ &=\frac {1}{2} \text {Subst}\left (\int \left (\frac {A}{a x^2}+\frac {-A b+a B}{a^2 x}+\frac {-a b B+A \left (b^2-a c\right )+(A b-a B) c x}{a^2 \left (a+b x+c x^2\right )}\right ) \, dx,x,x^2\right )\\ &=-\frac {A}{2 a x^2}-\frac {(A b-a B) \log (x)}{a^2}+\frac {\text {Subst}\left (\int \frac {-a b B+A \left (b^2-a c\right )+(A b-a B) c x}{a+b x+c x^2} \, dx,x,x^2\right )}{2 a^2}\\ &=-\frac {A}{2 a x^2}-\frac {(A b-a B) \log (x)}{a^2}+\frac {(A b-a B) \text {Subst}\left (\int \frac {b+2 c x}{a+b x+c x^2} \, dx,x,x^2\right )}{4 a^2}+\frac {\left (-a b B+A \left (b^2-2 a c\right )\right ) \text {Subst}\left (\int \frac {1}{a+b x+c x^2} \, dx,x,x^2\right )}{4 a^2}\\ &=-\frac {A}{2 a x^2}-\frac {(A b-a B) \log (x)}{a^2}+\frac {(A b-a B) \log \left (a+b x^2+c x^4\right )}{4 a^2}-\frac {\left (-a b B+A \left (b^2-2 a c\right )\right ) \text {Subst}\left (\int \frac {1}{b^2-4 a c-x^2} \, dx,x,b+2 c x^2\right )}{2 a^2}\\ &=-\frac {A}{2 a x^2}+\frac {\left (a b B-A \left (b^2-2 a c\right )\right ) \tanh ^{-1}\left (\frac {b+2 c x^2}{\sqrt {b^2-4 a c}}\right )}{2 a^2 \sqrt {b^2-4 a c}}-\frac {(A b-a B) \log (x)}{a^2}+\frac {(A b-a B) \log \left (a+b x^2+c x^4\right )}{4 a^2}\\ \end {align*}

________________________________________________________________________________________

Mathematica [A]
time = 0.15, size = 186, normalized size = 1.66 \begin {gather*} \frac {-\frac {2 a A}{x^2}+4 (-A b+a B) \log (x)+\frac {\left (-a B \left (b+\sqrt {b^2-4 a c}\right )+A \left (b^2-2 a c+b \sqrt {b^2-4 a c}\right )\right ) \log \left (b-\sqrt {b^2-4 a c}+2 c x^2\right )}{\sqrt {b^2-4 a c}}+\frac {\left (a B \left (b-\sqrt {b^2-4 a c}\right )+A \left (-b^2+2 a c+b \sqrt {b^2-4 a c}\right )\right ) \log \left (b+\sqrt {b^2-4 a c}+2 c x^2\right )}{\sqrt {b^2-4 a c}}}{4 a^2} \end {gather*}

Antiderivative was successfully verified.

[In]

Integrate[(A + B*x^2)/(x^3*(a + b*x^2 + c*x^4)),x]

[Out]

((-2*a*A)/x^2 + 4*(-(A*b) + a*B)*Log[x] + ((-(a*B*(b + Sqrt[b^2 - 4*a*c])) + A*(b^2 - 2*a*c + b*Sqrt[b^2 - 4*a
*c]))*Log[b - Sqrt[b^2 - 4*a*c] + 2*c*x^2])/Sqrt[b^2 - 4*a*c] + ((a*B*(b - Sqrt[b^2 - 4*a*c]) + A*(-b^2 + 2*a*
c + b*Sqrt[b^2 - 4*a*c]))*Log[b + Sqrt[b^2 - 4*a*c] + 2*c*x^2])/Sqrt[b^2 - 4*a*c])/(4*a^2)

________________________________________________________________________________________

Maple [A]
time = 0.05, size = 126, normalized size = 1.12

method	result	size
default	\(-\frac {\frac {\left (-b c A +a c B \right ) \ln \left (c \,x^{4}+b \,x^{2}+a \right )}{2 c}+\frac {2 \left (a c A -A \,b^{2}+a b B -\frac {\left (-b c A +a c B \right ) b}{2 c}\right ) \arctan \left (\frac {2 c \,x^{2}+b}{\sqrt {4 a c -b^{2}}}\right )}{\sqrt {4 a c -b^{2}}}}{2 a^{2}}-\frac {A}{2 a \,x^{2}}+\frac {\left (-A b +a B \right ) \ln \left (x \right )}{a^{2}}\)	\(126\)
risch	\(-\frac {A}{2 a \,x^{2}}-\frac {\ln \left (x \right ) A b}{a^{2}}+\frac {\ln \left (x \right ) B}{a}+\frac {\left (\munderset {\textit {\_R} =\RootOf \left (\left (4 a^{3} c -a^{2} b^{2}\right ) \textit {\_Z}^{2}+\left (-4 A a b c +A \,b^{3}+4 a^{2} c B -B a \,b^{2}\right ) \textit {\_Z} +A^{2} c^{2}-A B b c +B^{2} a c \right )}{\sum }\textit {\_R} \ln \left (\left (\left (-10 a^{3} c +3 a^{2} b^{2}\right ) \textit {\_R}^{2}+\left (4 A a b c -5 a^{2} c B \right ) \textit {\_R} -2 A^{2} c^{2}\right ) x^{2}+a^{3} b \,\textit {\_R}^{2}+\left (-a^{2} c A +2 A a \,b^{2}-2 B \,a^{2} b \right ) \textit {\_R} -2 A^{2} b c +2 A B a c \right )\right )}{2}\)	\(197\)

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((B*x^2+A)/x^3/(c*x^4+b*x^2+a),x,method=_RETURNVERBOSE)

[Out]

-1/2/a^2*(1/2*(-A*b*c+B*a*c)/c*ln(c*x^4+b*x^2+a)+2*(a*c*A-A*b^2+a*b*B-1/2*(-A*b*c+B*a*c)*b/c)/(4*a*c-b^2)^(1/2
)*arctan((2*c*x^2+b)/(4*a*c-b^2)^(1/2)))-1/2*A/a/x^2+1/a^2*(-A*b+B*a)*ln(x)

________________________________________________________________________________________

Maxima [F(-2)]
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {Exception raised: ValueError} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((B*x^2+A)/x^3/(c*x^4+b*x^2+a),x, algorithm="maxima")

[Out]

Exception raised: ValueError >> Computation failed since Maxima requested additional constraints; using the 'a
ssume' command before evaluation *may* help (example of legal syntax is 'assume(4*a*c-b^2>0)', see `assume?` f
or more deta

________________________________________________________________________________________

Fricas [A]
time = 0.48, size = 385, normalized size = 3.44 \begin {gather*} \left [\frac {{\left (B a b - A b^{2} + 2 \, A a c\right )} \sqrt {b^{2} - 4 \, a c} x^{2} \log \left (\frac {2 \, c^{2} x^{4} + 2 \, b c x^{2} + b^{2} - 2 \, a c + {\left (2 \, c x^{2} + b\right )} \sqrt {b^{2} - 4 \, a c}}{c x^{4} + b x^{2} + a}\right ) - 2 \, A a b^{2} + 8 \, A a^{2} c - {\left (B a b^{2} - A b^{3} - 4 \, {\left (B a^{2} - A a b\right )} c\right )} x^{2} \log \left (c x^{4} + b x^{2} + a\right ) + 4 \, {\left (B a b^{2} - A b^{3} - 4 \, {\left (B a^{2} - A a b\right )} c\right )} x^{2} \log \left (x\right )}{4 \, {\left (a^{2} b^{2} - 4 \, a^{3} c\right )} x^{2}}, \frac {2 \, {\left (B a b - A b^{2} + 2 \, A a c\right )} \sqrt {-b^{2} + 4 \, a c} x^{2} \arctan \left (-\frac {{\left (2 \, c x^{2} + b\right )} \sqrt {-b^{2} + 4 \, a c}}{b^{2} - 4 \, a c}\right ) - 2 \, A a b^{2} + 8 \, A a^{2} c - {\left (B a b^{2} - A b^{3} - 4 \, {\left (B a^{2} - A a b\right )} c\right )} x^{2} \log \left (c x^{4} + b x^{2} + a\right ) + 4 \, {\left (B a b^{2} - A b^{3} - 4 \, {\left (B a^{2} - A a b\right )} c\right )} x^{2} \log \left (x\right )}{4 \, {\left (a^{2} b^{2} - 4 \, a^{3} c\right )} x^{2}}\right ] \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((B*x^2+A)/x^3/(c*x^4+b*x^2+a),x, algorithm="fricas")

[Out]

[1/4*((B*a*b - A*b^2 + 2*A*a*c)*sqrt(b^2 - 4*a*c)*x^2*log((2*c^2*x^4 + 2*b*c*x^2 + b^2 - 2*a*c + (2*c*x^2 + b)
*sqrt(b^2 - 4*a*c))/(c*x^4 + b*x^2 + a)) - 2*A*a*b^2 + 8*A*a^2*c - (B*a*b^2 - A*b^3 - 4*(B*a^2 - A*a*b)*c)*x^2
*log(c*x^4 + b*x^2 + a) + 4*(B*a*b^2 - A*b^3 - 4*(B*a^2 - A*a*b)*c)*x^2*log(x))/((a^2*b^2 - 4*a^3*c)*x^2), 1/4
*(2*(B*a*b - A*b^2 + 2*A*a*c)*sqrt(-b^2 + 4*a*c)*x^2*arctan(-(2*c*x^2 + b)*sqrt(-b^2 + 4*a*c)/(b^2 - 4*a*c)) -
2*A*a*b^2 + 8*A*a^2*c - (B*a*b^2 - A*b^3 - 4*(B*a^2 - A*a*b)*c)*x^2*log(c*x^4 + b*x^2 + a) + 4*(B*a*b^2 - A*b
^3 - 4*(B*a^2 - A*a*b)*c)*x^2*log(x))/((a^2*b^2 - 4*a^3*c)*x^2)]

________________________________________________________________________________________

Sympy [F(-1)] Timed out
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {Timed out} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((B*x**2+A)/x**3/(c*x**4+b*x**2+a),x)

[Out]

Timed out

________________________________________________________________________________________

Giac [A]
time = 5.46, size = 124, normalized size = 1.11 \begin {gather*} -\frac {{\left (B a - A b\right )} \log \left (c x^{4} + b x^{2} + a\right )}{4 \, a^{2}} + \frac {{\left (B a - A b\right )} \log \left (x^{2}\right )}{2 \, a^{2}} - \frac {{\left (B a b - A b^{2} + 2 \, A a c\right )} \arctan \left (\frac {2 \, c x^{2} + b}{\sqrt {-b^{2} + 4 \, a c}}\right )}{2 \, \sqrt {-b^{2} + 4 \, a c} a^{2}} - \frac {B a x^{2} - A b x^{2} + A a}{2 \, a^{2} x^{2}} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((B*x^2+A)/x^3/(c*x^4+b*x^2+a),x, algorithm="giac")

[Out]

-1/4*(B*a - A*b)*log(c*x^4 + b*x^2 + a)/a^2 + 1/2*(B*a - A*b)*log(x^2)/a^2 - 1/2*(B*a*b - A*b^2 + 2*A*a*c)*arc
tan((2*c*x^2 + b)/sqrt(-b^2 + 4*a*c))/(sqrt(-b^2 + 4*a*c)*a^2) - 1/2*(B*a*x^2 - A*b*x^2 + A*a)/(a^2*x^2)

________________________________________________________________________________________

Mupad [B]
time = 4.85, size = 2500, normalized size = 22.32 \begin {gather*} \text {Too large to display} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((A + B*x^2)/(x^3*(a + b*x^2 + c*x^4)),x)

[Out]

- A/(2*a*x^2) - (log(x)*(A*b - B*a))/a^2 - (log(((A^3*c^5*x^2)/a^3 - (((((4*b*c^2*(A*a*c - A*b^2 + B*a*b))/a -
(2*c^3*x^2*(A*b^2 + 10*A*a*c - 5*B*a*b))/a + (b*c^2*(a*b + 3*b^2*x^2 - 10*a*c*x^2)*(B*a - A*b + a^2*(-(2*A*a*
c - A*b^2 + B*a*b)^2/(a^4*(4*a*c - b^2)))^(1/2)))/a^2)*(B*a - A*b + a^2*(-(2*A*a*c - A*b^2 + B*a*b)^2/(a^4*(4*
a*c - b^2)))^(1/2)))/(4*a^2) + (A*c^3*(A*a*c - 4*A*b^2 + 4*B*a*b))/a^2 - (A*c^4*x^2*(6*A*b - 5*B*a))/a^2)*(B*a
- A*b + a^2*(-(2*A*a*c - A*b^2 + B*a*b)^2/(a^4*(4*a*c - b^2)))^(1/2)))/(4*a^2) + (A^2*c^4*(A*b - B*a))/a^3)*(
(((((2*c^3*x^2*(A*b^2 + 10*A*a*c - 5*B*a*b))/a - (4*b*c^2*(A*a*c - A*b^2 + B*a*b))/a + (b*c^2*(a*b + 3*b^2*x^2
- 10*a*c*x^2)*(A*b - B*a + a^2*(-(2*A*a*c - A*b^2 + B*a*b)^2/(a^4*(4*a*c - b^2)))^(1/2)))/a^2)*(A*b - B*a + a
^2*(-(2*A*a*c - A*b^2 + B*a*b)^2/(a^4*(4*a*c - b^2)))^(1/2)))/(4*a^2) + (A*c^3*(A*a*c - 4*A*b^2 + 4*B*a*b))/a^
2 - (A*c^4*x^2*(6*A*b - 5*B*a))/a^2)*(A*b - B*a + a^2*(-(2*A*a*c - A*b^2 + B*a*b)^2/(a^4*(4*a*c - b^2)))^(1/2)
))/(4*a^2) + (A^3*c^5*x^2)/a^3 + (A^2*c^4*(A*b - B*a))/a^3))*(2*A*b^3 - 2*B*a*b^2 + 8*B*a^2*c - 8*A*a*b*c))/(2
*(16*a^3*c - 4*a^2*b^2)) - (atan((16*a^6*x^2*((((((5*A*B*a^2*c^4 - 6*A^2*a*b*c^4)/a^3 - (((20*A*a^3*c^4 - 10*B
*a^3*b*c^3 + 2*A*a^2*b^2*c^3)/a^3 + ((40*a^4*b*c^3 - 12*a^3*b^3*c^2)*(2*A*b^3 - 2*B*a*b^2 + 8*B*a^2*c - 8*A*a*
b*c))/(2*a^3*(16*a^3*c - 4*a^2*b^2)))*(2*A*b^3 - 2*B*a*b^2 + 8*B*a^2*c - 8*A*a*b*c))/(2*(16*a^3*c - 4*a^2*b^2)
))*(2*A*b^3 - 2*B*a*b^2 + 8*B*a^2*c - 8*A*a*b*c))/(2*(16*a^3*c - 4*a^2*b^2)) - (A^3*c^5)/a^3 + (((((20*A*a^3*c
^4 - 10*B*a^3*b*c^3 + 2*A*a^2*b^2*c^3)/a^3 + ((40*a^4*b*c^3 - 12*a^3*b^3*c^2)*(2*A*b^3 - 2*B*a*b^2 + 8*B*a^2*c
- 8*A*a*b*c))/(2*a^3*(16*a^3*c - 4*a^2*b^2)))*(2*A*a*c - A*b^2 + B*a*b))/(4*a^2*(4*a*c - b^2)^(1/2)) + ((40*a
^4*b*c^3 - 12*a^3*b^3*c^2)*(2*A*a*c - A*b^2 + B*a*b)*(2*A*b^3 - 2*B*a*b^2 + 8*B*a^2*c - 8*A*a*b*c))/(8*a^5*(4*
a*c - b^2)^(1/2)*(16*a^3*c - 4*a^2*b^2)))*(2*A*a*c - A*b^2 + B*a*b))/(4*a^2*(4*a*c - b^2)^(1/2)) + ((40*a^4*b*
c^3 - 12*a^3*b^3*c^2)*(2*A*a*c - A*b^2 + B*a*b)^2*(2*A*b^3 - 2*B*a*b^2 + 8*B*a^2*c - 8*A*a*b*c))/(32*a^7*(4*a*
c - b^2)*(16*a^3*c - 4*a^2*b^2)))*(3*A*b^4 + A*a^2*c^2 - 3*B*a*b^3 - 9*A*a*b^2*c + 8*B*a^2*b*c))/(8*a^3*c^2*(2
5*B^2*a^3*c - 6*A^2*b^4 + A^2*a^2*c^2 - 6*B^2*a^2*b^2 + 12*A*B*a*b^3 + 24*A^2*a*b^2*c - 49*A*B*a^2*b*c)) + (((
((5*A*B*a^2*c^4 - 6*A^2*a*b*c^4)/a^3 - (((20*A*a^3*c^4 - 10*B*a^3*b*c^3 + 2*A*a^2*b^2*c^3)/a^3 + ((40*a^4*b*c^
3 - 12*a^3*b^3*c^2)*(2*A*b^3 - 2*B*a*b^2 + 8*B*a^2*c - 8*A*a*b*c))/(2*a^3*(16*a^3*c - 4*a^2*b^2)))*(2*A*b^3 -
2*B*a*b^2 + 8*B*a^2*c - 8*A*a*b*c))/(2*(16*a^3*c - 4*a^2*b^2)))*(2*A*a*c - A*b^2 + B*a*b))/(4*a^2*(4*a*c - b^2
)^(1/2)) - (((((20*A*a^3*c^4 - 10*B*a^3*b*c^3 + 2*A*a^2*b^2*c^3)/a^3 + ((40*a^4*b*c^3 - 12*a^3*b^3*c^2)*(2*A*b
^3 - 2*B*a*b^2 + 8*B*a^2*c - 8*A*a*b*c))/(2*a^3*(16*a^3*c - 4*a^2*b^2)))*(2*A*a*c - A*b^2 + B*a*b))/(4*a^2*(4*
a*c - b^2)^(1/2)) + ((40*a^4*b*c^3 - 12*a^3*b^3*c^2)*(2*A*a*c - A*b^2 + B*a*b)*(2*A*b^3 - 2*B*a*b^2 + 8*B*a^2*
c - 8*A*a*b*c))/(8*a^5*(4*a*c - b^2)^(1/2)*(16*a^3*c - 4*a^2*b^2)))*(2*A*b^3 - 2*B*a*b^2 + 8*B*a^2*c - 8*A*a*b
*c))/(2*(16*a^3*c - 4*a^2*b^2)) + ((40*a^4*b*c^3 - 12*a^3*b^3*c^2)*(2*A*a*c - A*b^2 + B*a*b)^3)/(64*a^9*(4*a*c
- b^2)^(3/2)))*(6*A*b^5 - 20*B*a^3*c^2 - 6*B*a*b^4 - 30*A*a*b^3*c + 26*A*a^2*b*c^2 + 28*B*a^2*b^2*c))/(16*a^3
*c^2*(4*a*c - b^2)^(1/2)*(25*B^2*a^3*c - 6*A^2*b^4 + A^2*a^2*c^2 - 6*B^2*a^2*b^2 + 12*A*B*a*b^3 + 24*A^2*a*b^2
*c - 49*A*B*a^2*b*c)))*(4*a*c - b^2)^(3/2))/(4*A^2*a^2*c^4 + A^2*b^4*c^2 + B^2*a^2*b^2*c^2 - 4*A^2*a*b^2*c^3 -
2*A*B*a*b^3*c^2 + 4*A*B*a^2*b*c^3) + (a^3*(4*a*c - b^2)*((((((4*A*a^3*b*c^3 - 4*A*a^2*b^3*c^2 + 4*B*a^3*b^2*c
^2)/a^3 + (2*a*b^2*c^2*(2*A*b^3 - 2*B*a*b^2 + 8*B*a^2*c - 8*A*a*b*c))/(16*a^3*c - 4*a^2*b^2))*(2*A*a*c - A*b^2
+ B*a*b))/(4*a^2*(4*a*c - b^2)^(1/2)) + (b^2*c^2*(2*A*a*c - A*b^2 + B*a*b)*(2*A*b^3 - 2*B*a*b^2 + 8*B*a^2*c -
8*A*a*b*c))/(2*a*(4*a*c - b^2)^(1/2)*(16*a^3*c - 4*a^2*b^2)))*(2*A*b^3 - 2*B*a*b^2 + 8*B*a^2*c - 8*A*a*b*c))/
(2*(16*a^3*c - 4*a^2*b^2)) + (((A^2*a^2*c^4 - 4*A^2*a*b^2*c^3 + 4*A*B*a^2*b*c^3)/a^3 + (((4*A*a^3*b*c^3 - 4*A*
a^2*b^3*c^2 + 4*B*a^3*b^2*c^2)/a^3 + (2*a*b^2*c^2*(2*A*b^3 - 2*B*a*b^2 + 8*B*a^2*c - 8*A*a*b*c))/(16*a^3*c - 4
*a^2*b^2))*(2*A*b^3 - 2*B*a*b^2 + 8*B*a^2*c - 8*A*a*b*c))/(2*(16*a^3*c - 4*a^2*b^2)))*(2*A*a*c - A*b^2 + B*a*b
))/(4*a^2*(4*a*c - b^2)^(1/2)) - (b^2*c^2*(2*A*a*c - A*b^2 + B*a*b)^3)/(16*a^5*(4*a*c - b^2)^(3/2)))*(6*A*b^5
- 20*B*a^3*c^2 - 6*B*a*b^4 - 30*A*a*b^3*c + 26*A*a^2*b*c^2 + 28*B*a^2*b^2*c))/(c^2*(4*A^2*a^2*c^4 + A^2*b^4*c^
2 + B^2*a^2*b^2*c^2 - 4*A^2*a*b^2*c^3 - 2*A*B*a*b^3*c^2 + 4*A*B*a^2*b*c^3)*(25*B^2*a^3*c - 6*A^2*b^4 + A^2*a^2
*c^2 - 6*B^2*a^2*b^2 + 12*A*B*a*b^3 + 24*A^2*a*b^2*c - 49*A*B*a^2*b*c)) - (2*a^3*(4*a*c - b^2)^(3/2)*((A^3*b*c
^4 - A^2*B*a*c^4)/a^3 - (((A^2*a^2*c^4 - 4*A^2*a*b^2*c^3 + 4*A*B*a^2*b*c^3)/a^3 + (((4*A*a^3*b*c^3 - 4*A*a^2*b
^3*c^2 + 4*B*a^3*b^2*c^2)/a^3 + (2*a*b^2*c^2*(2*A*b^3 - 2*B*a*b^2 + 8*B*a^2*c - 8*A*a*b*c))/(16*a^3*c - 4*a^2*
b^2))*(2*A*b^3 - 2*B*a*b^2 + 8*B*a^2*c - 8*A*a*...

________________________________________________________________________________________